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ABSTRACT 

The recent discovery of day-long gamma-ray flares in the Crab Nebula, pre- 
sumed to be synchrotron emission by PeV (10 15 eV) electrons in milligauss mag- 
netic fields, presents a strong challenge to particle acceleration models. The ob- 
served photon energies exceed the upper limit (~100 MeV) obtained by balancing 
the acceleration rate and synchrotron radiation losses under standard conditions 
where the electric field is smaller than the magnetic field. We argue that a linear 
electric accelerator, operating at magnetic reconnection sites, is able to circum- 
vent this difficulty. Sufficiently energetic electrons have gyroradii so large that 
their motion is insensitive to small-scale turbulent structures in the reconnection 
layer and is controlled only by large-scale fields. We show that such particles 
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are guided into the reconnection layer by the reversing magnetic field as they are 
accelerated by the reconnection electric field. As these electrons become confined 
within the current sheet, they experience a decreasing perpendicular magnetic 
field that may drop below the accelerating electric field. This enables them to 
reach higher energies before suffering radiation losses and hence to emit syn- 
chrotron radiation in excess of the 100 MeV limit, providing a natural resolution 
to the Crab gamma-ray flare paradox. 

Subject headings: Acceleration of particles — magnetic reconnection — radia- 
tion mechanisms: non-thermal — pulsars: individual (Crab) — ISM: individual 
objects (Crab nebula) — gamma rays: stars 



1. Introduction 



Rece n tly discovered high -energy (> 100 MeV) gamma-ray flares in the Crab Nebula (ITavani et al. 



2011 



tion f Lyutikov 



2010 



Abdo et al.l |2011|) challenge theoretical mod e ls of particle accelerat ion and radia 



Komissarov fc Lyutikovl |2010| ; iBednarek fc Idea 120101 ) . These flares 



are believed to come from synch rotron radiation in t he nebula, i.e., outside t he pulsar wind 
termination shock, at > 10 16 cm ( iRees &: Gunnlll974l ; lKennel fc Coronitilll984l ). and not from 
the vicinity of the pulsar itself. Their most puzzling aspect is the very high photon energies, 
reaching a few hundred MeV. This exceeds the accepted upper limit (radiation-reaction limit) 
for photons produced via the synchrotron mechanism, e syncc ~ (9/4)m e c 2 /ttf s ~ 160 MeV 
(where = e 2 /hc is the fine structure constant), obtained by balancing the acceler- 
ating electric force on a particle eE with the synchrotron radiation reaction drag force, 
/ rad ~ Psync/c = 2cr< r 7 2 -E>i/87r (where 7 is the Lorentz factor of the particle, <jt = 87re 4 /3mgC 4 
is the Thomson cross-section and "_L" means perpendicular to the particle's motion), in com- 
bination with the requirement E < B±. The latter condition is required in most established 
particle acceleration mechanisms, e.g., diffusive shock acceleration; thus, these mechanisms 
cannot ac celerate particles to energies high enough to produce the observed synchrotron 
radiation (IGuilbert et al.lll983l ; Ide Jager et al.lll996l ). 



A possible resolution of this para dox is to invoke relativistic Doppler boosting of the emitting 
regio n towards the observer (e.g. 



Lvutikov 


2010 


Komissarov & Lvutikov 


2010 


Bednarek & Idee 



2010h . However, although the required bulk Doppler factors (~ 3) are not ruled out, 
there is lit tle evidence for su ch high speeds from proper motions in the inner regions of 
the Crab (IHester et al.l 120021 ) and theoretical arguments suggest that typical bulk speeds 
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dow nstream of the shock are at be st mildly relativistic, unless the shock is strongly oblique 
(see 



Komissarov fc Lyutikovll2010l ) 



Here we focus on the alternative possibility that electrons are indeed accelerated to very 
high (PeV) energies in an extended region where E > B±. This is impossible in ideal mag- 
netohydrodynamics (MHD), E = — v x B/c, where v < c is the bulk flow velocity, implying 
E < B± (assuming that the particle moves along E). However, ideal MHD breaks down at 
magnetic reconnection sites, where E > B± is satisfied naturally in a thin layer. Sufficiently 
energetic particles tend to be drawn into the layer, wh ere they be come trapped and may 
be accelerated to extreme energies, as was suggested by iKirki ( 120041 ) . By explicitly demon- 
strating this focusing effect in the ultrarelativistic limit, in this paper we show that a linear 
accelerator, utilizing the large-scale electric field associated with a reconnection process, can 
explain gamma-ray flares in the Crab. We describe the mechanism in § 2 and its application 
to the Crab flares in § 3. 



2. Extreme Particle Acceleration and Radiation in Reconnecting Current 

Sheets 



Reconnection is generally recognized as an important mechanism of non-thermal particle 
accel eration, including in relativistic pair plasmas thought to exist in the Crab Nebula 



(e.g., iZenitani fc Hoshinoll200ll 120081 ) . Acceleration mechanisms at moderate energies may 
be co mplex, e.g., involvi ng small-scale turbulent structures such as multiple magnetic is- 
lands (IDrake et al.l 120061 1 . However, the most energetic particles (e.g., the PeV particles 
considered here) are special: their relativistic gyroradii, ri = ^p c ~ 1.7 x 10 15 cm £^379 
(where B^ 3 = B/10~ 3 G, 79 = 7/IO 9 , and p c = c/co c = m e c 2 /eB ~ 17km Bz\) are large, 
comparable to the global flare region size I (< 1 light day ~ 3 x 10 15 cm based on observed 
flare duration). Therefore, small-scale turbulence is washed out for these particles and the 
complicated mechanisms operating at moderate energies are not important for them. Their 
motion is sensitive only to the large-scale electromagnetic field structure (Figured]), which 
consists of the reconnecting magnetic field B x (y), reversing across a current layer of some 
thickness 5; the nearly- uniform reconnection electric field E z \ and perhaps a guide mag- 
netic field B z , also roughly uniform. A finite B z may in fact be required for accelerat ion 
of moderate-energy particles in relativistic pair reconnection ( iZenitani &: Hoshinol 120081 ). a 
pre-condition for injecting seed particles for the highest-energy acceleration considered in 
this paper. We neglect the reconnected magnetic field B y . (Note that the Hall magnetic and 
electric fields are absent in pair plasmas.) 
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In this electromagnetic field, ultra-rel ativistic particles move mainly in the ^ - direction i n 
a relativistic analog of Speiser orbits (jSpeiserl 11965c IZenitani fc Hoshinoll2001t iKirkl 120041 ) . 
They are accelerated by the electric field E z in the ^-direction while being confined to 
the layer midplane (y = 0) by the reversing reconnecting magnetic field B x . In the re- 
direction, they may be confined by the guide magnetic field. Importantly, as the particles get 
accelerated, their orbits not only stretch in the ^-direction, but also shrink in the ^/-direction; 
both the Speiser meandering width y max (the maximum deviation from the midplane) and 
the midplane-crossing angle 8q decrease as their energy increases. Thus, the population of 
energetic particles, even if isotropic initially, gets focused into a tight beam along the current 
sheet. 



We explain this focusing mechanism with the following analytical model. For simplicity, 
assume B z = and consider an electron whose trajectory lies entirely in the (yz) plane. We 
choose the starting point at the origin (0,0,0), so that at z = the particle is just crossing 
the midplane y = at some angle 9q < 7r /2 (see Figured]). We also assume that the electron 
is already pre-accelerated to a very high Lorentz factor 70 3> 1 and that 8q is not too small, 
so that both its gyroradius and y max are much larger than the current layer thickness 5; then, 
most of the trajectory lies outside the layer, in the two upstream regions with \B x (y)\ ps B = 
const. The electric field is uniform, E z = — E = — f3 rcc B ~ — 0.3 V/cm/3 rec B _ 3 , where 
the dimensi qnless reconnection rate /3 ren < 1 is ~ 0.1, as indicated by relativistic particle 
simulations (IZenitani fc Hoshinoll2008l ). Finally, we take the initial particle energy to be well 
below the radiation reaction limit. 



To determine the evolution of the orbit's key parameters, first consider one half-cycle segment 
of the orbit, z G [0, zi], where Z\ is the next midplane-crossing. The motion is determined 
by two parameters: 70 and the midplane-crossing angle 8q. We wish to calculate the slight 
changes in 7 and \8\ over this segment, i.e., the values 7(^1) and from this we will 

derive the secular evolution of y mSLX . 

The 7-factor follows from energy conservation: e(z) = e + eE z, i.e., 

l{z) = 70 + PmcZ , (1) 

where z = zj p c . The particle's trajectory is obtained from the y component of the relativistic 
equation of motion: d( , jVy)/dt = —u c v z =3- d{ r yv y )/dz = —u c = —eB /m e c, integrating which 
we get 

7(2) p y (z) = 'joPyo - z ~ 70 sin 9 - z , (2) 

where f3 y o = (3 sin^o ~ sin^o since (3 = v/c ~ 1. This gives us the trajectory's apex z a 
(where v y = and y = y max ): z a = 70 sin6» . 
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Next, we assume that \9\ < 9q <C 1, so that dy/dz = f3 y J (3 Z ~ (3 y , where y = y/p c - Then, 
integrating Equation ([2]), we obtain an explicit expression for the particle trajectory: y(z) = 
— zj f3 TCC + (7o//^rcc) (1 + /^rccSin^o) ln(l + /3 rec ^/7o)- Anticipating the fractional increase of 
the particle's energy over z 1 to be small, /3 rcc ^i/7o ~ $o ^ 1 ; we expand y{z) and find 
y max = y{z a ) ~ 7o#o/2 and the segment's length z\ (the distance to the next midplane 
crossing, defined by y{z\) = 0): z\ ~ 270^0 + (2/3) 7o/3 rC c#o- This allows us to estimate 
the changes in the orbit parameters, 7 and |#o|> from one midplane crossing to the next: 
5 1 = 7 (*l) -7o = &ec% ^ 2/3 rcc7( A > 0, and 5\9 \ = \0{z x )\ - 6 ~ - (4/3) /3 rcc # 2 < 0. 
Thus, as 7 increases, |#o| decreases, i.e., the trajectory becomes increasingly aligned with the 
accelerating electric field. 

The secular evolution of 7 and \8 \ over many current-sheet crossing cycles (z ^> Z\) follows 
from d\9o\/d^ ~ 8\9 G \/8^ ~ —(2/3) |^o|/t? integrating which we get \0 Q \ ~ 7~ 2 / 3 ~ z~ 2 / 3 . 
Similarly, we find z\ ~ 7 1 / 3 ~ z 1 ^ 3 and ?/ max ~ 7 -1 / 3 ~ 2" 1 / 3 . Thus, the highest-energy 
particles are focused into a narrow beam confined closer and closer to the midplane! (This 
shrinkage of the trajectory can be interpreted as a result of the E x B drift of the particle's 
virtual guiding center away from the midplane, separately in each half-cycle segment. It can 
also be interpreted dbS ci result of the conservation of the cidicibcitic invciricint Jy — J Pydy ~ 
7m e c6» y max oc 7 2 0g = const.) 

Two factors limit the particle energy. First, the finite current-sheet length, I = 10 16 /i6 cm ; 
limits the energy to 

e m ax = 7max m e c 2 = eE l = e(3 rec B 1 ~ 3 PeV f3 TCC B _ 3 1 16 , (3) 

corresponding to 7 max = 6 x 10 9 f3 rcc B ^ 3 / 16 and 

esyncmax = (3/2) 7max = (3/2) (I / p c ) 2 hw c ~ 600 MeV /3 2 cc £ 3 _ 3 / 2 6 . (4) 



Second, radiation reaction may cause the energy to saturate at a lower value. If the initial 
injection values y max , in j and 7 inj satisfy y max>ini /5 > (7rad,*/7inj) 1/3 (where 7 rad) * is defined in 
Equation ([5]) below), then the radiation reaction-limited regime is reached while most of the 
orbit is outside the layer, where \B\ = B . The cycle- averaged radiation reaction force / ra d 
balances the electric force when 

7 = 7rad,* = J^-Prec ^ 3 X 1 9 fig B~^ . (5) 

By comparing 7 max and 7 ra d,*, we see that the radiation reaction limit is reached before 7 max 
only if the layer is sufficiently long, Z 16 > 0.5f3^ 2 Bq 3 / 2 . Then, the particle moves in a 
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sinusoidal-like orbit with 7 ~ 7 ra d,o = const and efficiently converts all the energy it gains 
from the electric field into synchrotron photons with a characteristic energy 

e sy nc* = (3/2) 7r 2 adi , hw c = (9/4) ft ec c^ 1 m e c 2 = 160 MeV (3 rec , (6) 

somewhat too low to explain the Crab flares. 

However, despite the saturation of 7, the particle orbit continues to shrink towards the mid- 
plane with 6*o ~ exp(— /3 rcc 2/3p c7ra d,*). If the layer is long enough, then ?/ max eventually 
shrinks below 5, reducing radiative losses and thus making extremely high particle energies 
possible. Indeed, when the entire trajectory is contained deep within the layer, the maxi- 
mum perpendicular magnetic field sampled by the particle, -B max = B x (y max ) ~ B y miix /S, 
becomes smaller than Bq, thus reducing / ra d for a fixed 7, whereas the accelerating force 
eEo = e/3 rec -Bo remains unchanged. Correspondingly, the radiation reaction limit 7ra d in- 
creases in compensation to maintain the balance between the cycle-averaged / ra d and the 
constant electric acceleration: 7 rad « 7ra d,* ^/Umax., an d may rise well above 7 ra d,* as y max 
shrinks. The corresponding synchrotron photon energy, e sync ~ e sy nc,* Vj/maxj can eas ily 
exceed e synC) *, and thus can explain the > 100 MeV photon energies observed in the Crab 
flares. 

We can derive a relationship between 9q and 7 ra d describing their joint evolution in this 
regime. As a rough estimate, we can regard -B max as playing the same role as Bq played 
for a particle for which most of the trajectory lied outside the layer. Then we can es- 
timate the meandering width as y max ~ 0.5 7rad 9% p c B /B max ~ 0.5 7rad 9% p c 5/y max , i.e., 
y max ~ 9 ("frndPcd/Z) 1 / 2 - Combining this with the expression for the radiation-reaction 
limited Lorentz factor 7 rad ~ 7ra d,*^/2/m a x where 7ra d,* is given by Equation (jSJ), we find 
^o^7 r ad,*(25/Pc) 1/2 7~f. 

Alternatively, if the initial y max> i n j an d 7 i n j are small enough, then the orbit's meander- 
ing width |/ max quickly collapses below 5, so that the radiation reaction does not become 
important before the particle reaches the end of the layer. In this case, the particle is sim- 
ply accelerated by the electric field almost up to the theoretical maximum limit 7max = 
6 x 10 9 j3 mc -80,-3 ^16; which may exceed the conventional radiation-reaction limit given by 
Equation (j3J). It will then radiate all its energy in one short (compared to its relativistic 
cyclotron period 2t[^uj~ 1 ) and powerful burst when it finally escapes the reconnection layer 
into an adjacent region with a finite magnetic field {B ~ Bq), i.e., when it transitions from 
a big linear accelerator to a powerful radiator (beam dump). If this transition is sharp, then 
the characteristic synchrotron photon energy will be that given by Equation (j4j). 
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As ultrarelativistic particles moving on relativistic Speiser orbits cross the layer midplane, 
they effectively see a rapidly reversing magnetic field, with a characteristic reversal scale 
As ~ Z\ much smaller than their gyroradius — ip r .. This raises the question of whether 
the resulting radiation could in fact be jitter radiation ( lMedvedevll2000l ). However, the jitter 
parameter fitter = ^b/Pc is much larger than 1 in our case; for example, for particles with 
Umax 3> 5, we have \b ~ z\ ~ 2 r y 9 p c 3> p c . Therefore, the radiation emitted by these 
particles is true synchrotron radiation. 



Our analytical model is fully confirmed and generalized by our numerical calculation of 3D 
relativistic particle orbits (Figures |2]|4|) for arbitr ary initial pa rticle parameters, utilizing 
an explicit 8th-order Runge-Kutta-Verner method (jVernerlll978l ). We include the radiation 
reaction force, a guide magnetic field B z , and a realistic Harris profile of the reconnecting 
magnetic field, B x (y) = _B tanh {y/S). Physically, a lower limit on the current layer thick- 
ness S is set by the collisionless skin-depth or the relativistic gyroradius of the bulk electrons: 
<5 > 7buikPc = 1.7 x 10 12 cm7 bulk>6J Bl3, where 7 bu ik,6 = 7buik/10 6 . In reality, however, the 
layer may be broadened by turbulence driven by secondary instabilities. 



Our numerical study confirms the focusing effect and the above analytical scaling relation- 
ships (Cerutti et al. 2011, in preparation). Figure [2] shows an electron's orbit in the (yz)- 
plane with no guide field and the evolution of its Lorentz factor and its midplane-crossing 
angle along the ^-direction. A non-zero guide field adds a circular motion to the particle's 
orbit in the (xy)-pleaie (see Figure [3]), but the overall evolution of 7 and 9q is almost un- 
changed even for B z = Bq. Figure H] gives the characteristic energy e sync of synchrotron 
photons radiated by the electrons leaving the layer after four days of acceleration, as a func- 
tion of the inital parameters 7i n j and 6 in j in a 5 mG magnetic field. We find that > 100 MeV 
photons can be emitted for 7 in j < 10 7 regardless of the initial angle 9-^. The particles are 
confined within a tight cone of semi-aperture angle of only a few degrees. 



3. Application to Crab Gamma- Ray Flares 



In order to explain the short flare timescales (rg < 1 day, see lBalbo et al.ll201ll ). the flaring 
region must be compact (/ < 10 16 cm), suggesting that the flares originate just outside the 
pulsar wind termination shock. There are two places where reconnection current sheets may 
form in the inner Crab Nebula. First, the equator ial plane contains a reversing toroidal field 



flBegelmanl Il999l ; iKomissarov fc Lyubarsky 



20041 ). Second, we expect a cylindrical z-pinch 



along the rotational axis, which is unstable to kink-like MHD instabilities the nonlinear 
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development of which may lead to episodes of reconnection-driven dissipation (similar to 
qua si-periodic sawto o th crashes in tokam aks), perhaps powering X-ray emission in the axial 
jet (jBegelmanl Il998l ; iMizuno et al.ll201lf ). Since the PeV particles are focused along the 
current sheet, we can see the flare only when the reconnection electric field points toward 
us. While neither the equatorial plane nor the rotational axis lie close to our line of sight, 
instabilities may swing the current sheets through wide angles, allowing us to observe the 
radiation intermittently. 



Given the compact scale of the flaring region, the voltage drop EqI required to explain the 
energies of radiating electrons implies Bq > several mG. A similar estimate is obtained 
by requiring ta to exceed the synchrotron cooling time for > 100 MeV-emitting particles 
(synchrotron cooling can indeed be shown to dominate over adiabatic cooling), yielding 
-Bo -3 r fl,5 > 3; where r flj5 = r fl /10 5 s. These values are substantially higher th an the typical 



fields in the inner Crab Nebul a usually inferred from dynamical arguments (IRees fc Gunn 
1974 ; iKennel &: Coronit ill 19841 ) . But this discrepancy need not be problematic. First, if the 
axial z-pinch is responsible for the flares, then the field might be amplified by the pinch 
effect. Furthermore, the traditional estimates for the field strength in the nebula may be 
incorrect. Indeed, as is well-established, if ideal MHD holds everywhere outside the termi- 
nation shock, then the Poynting flux in the wind must be < 1% o f the kinetic energy flux 
( IRees fc Gunnlll974l ; IKennel fc Coronitilll984l ; iBegelman &: Lilll992[ ). This presents a prob- 
lem, since ultrarelativis tic, magnetically-drive n winds do not convert most of their magnetic 
energy to kinetic form ( IBegelman Sz Lilll994f ). Efficient magnetic reconnection downstream 
of the shock can reso lve this problem by dissipating m agnetic energy pumped into the nebula 
by the pulsar wind (jBegelmanl Il998l ; iLyutikovl 120101 ) . If this is the case, the flares provide 
an observational clue reconciling the physics of the pulsar wind with the properties of the 
nebula. 



The strong emission anisotropy in our model, with two oppositely directed beams (one 
produced by electrons and the other by positrons), can alleviate potential problems with the 
energetic efficiency of the flares. If the rad iation were isotropic, the total energ y of > 100 
MeV photons would be about 4 x 10 40 erg jTavani et aDhoilt kbdo et al.lboilh . Since the 



flare duration limits the emitting region si ze to about 1 light-day ~ 3 x 10 cm ( iTavani et al. 



201 ll ; lAbdo et al.ll201ll ; iBalbo et all 120111 ). and hence the volume to about 3 x 10 46 cm 3 , the 
required energy density deposited into PeV electrons (emitting 100 MeV photons) would 
have to be about 10~ 6 erg/cm 3 . If the energy source for the flare is magnetic (as it is in 
reconnection), then the required minimum field strength would be B min ~ 5mG K 1 / 2 , where 
K~ x < 1 is the fraction of the dissipated energy going to the PeV particles. Thus, even for 
a magnetic field of several mG, about 10 times higher than the standard value for the field 
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in the nebula, isotropic emission would require K ~ 1, i.e., a large fraction of the energy to 
go to PeV particles, which is problematic. However, having the emission beamed into two 
small solid angle cones alleviates this difficulty. 

Reconnection-driven particle acceleration may be a persistent mode of dissipation in the 
Crab. Flaring may occur all the time, but we see powerful gamma-ray flares only rarely, 
when the emission is beamed towards us. The strongest flares are observed only about twice 
a year, which is consistent with the typical single flare duration of a couple of days if the 
beam's solid angle is ~ 0.1 sr. Although we nominally associate the flare duration with 
the size of the acceleration region and the synchrotron cooling time for the most energetic 
electrons, it is also possible that the flare's intrinsic timescale is actually longer but that 
the current sheet and hence the beam wiggle around. Then, the observed time scale just 
corresponds to the time for the beam to cross our line of sight and one should expect a nearly 
symmetric light-curve profile. 

This mechanism may also be applicable to other astrophysica l systems, e.g., pulsar striped 



winds and blazar jets, as was first pointed out by iKirkl (120041 ) 
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Fig. 1. — Sketch of a relativistic Speiser orbit of a particle in a reconnection layer. 
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Fig. 2. — Numerically calculated orbit of a relativistic electron in reconnection layer of width 
8 = 7buikPc ~ 3.4 x 10 11 cm (y = ±5 is shown by red dotted lines), B = 5 mG, B z = 0, and 
Prec = 0.1. The particle is initially injected at the origin with 7 in j = 3 x 10 6 , # in j = 90°. The 
inserts describe the evolution of the particle's Lorentz factor 7 and of the midplane-crossing 
angle 9 . The vertical dotted line shows the distance z where the orbit becomes contained 
in the current layer. 
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Fig. 3. — Three-dimensional orbit of a relativistic electron in the reconnection layer with a 
non-zero guide field B z = B = 5 mG. The other parameters are 7i n j = 3 x 10 6 , 9 [n j = 90°, 
7buik = 10 6 and (3 rec = 0.1. The gray lines show the orbit projected onto the (xy) and (xz) 
planes. 
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Fig. 4. — Characteristic synchrotron photon energy produced by an electron at the end of 
the reconnection layer (7 = 4 light days), as a function of the particle's two initial parameters 
7; n j and 6^, with B = 5 mG and 5 = 3.4 x 10 11 cm. 



